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Abstract 

As voltage decreases d.c. condctivity of a Tomonaga-Luttinger liquid wire 
collapses to a small value determined by the length of the wire and its contacts 
with the leads. In condition that voltage drop (V) mostly occurs across a 
tunnel barrier inside the wire the tunneling density of states and, hence, 
the differential conductivity are shown to exhibit an interference structure 
resulted from the transition of the Luttinger liquid quasiparticles into free 
electrons at the exits from the wire. The finite length correction to the scale- 
invariant V 2 / 9 ~ 2 dependence of the conductivity oscillates as a function of 

voltage with periodicities related to both rigth and left traversal times. 
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Recently quantum transport in Tomonaga-Luttinger (TL) liquids has attracted a great 
deal of interest |I] as it was suggested to model both transport through a ID constriction 
|2|j3|] and the edge state transport in the fractional quantum Hall regime |3Q. The problem 
most often discussed in this context was a suppression of the TL transport by a point 
scatterer |3|,[|]. Finally, its exact solution has been constructed || which, however, could 
not be addressed to the ID transport to full extent because of the importance of the finite 
length effect in this case To handle the latter effect on the transport between two leads 
through a ID wire the model of the ihomogeneous TL liquid (ITL) has been invented PHTT 



It has allowed in particular to predict an interference structure in the density of states of 
the clean wire at low energy [p~2]j . 

The aim of my present work is to examine the finite length effect on the low voltage 
conductivity suppression in the ITL model of the one channel wire with a point scatterer 
inside. Temperature dependence of the conductance in this model was considered in |L3| . 



In spite of simplicity the model could be relevant to experiment due to its local stability 
against the influence of other impurities. Indeed, suppression of the current by an impurity 
is strengthened in (tiEp) 9-1 times in the TL model due to interaction. ( Here ti is the 
time of travelling from the impurity to the closest end of the wire, Ep is the Fermi energy 
in the wire and g < 1 is the constant of forward scattering.) It means that the effect of the 
impurities located closer to the middle of the wire is more important. On the other hand, 
among a few closely located impurities only one could bring about the strong suppression 
since in that case it cuts the wire and the renormalization of the scattering strength of the 
other impurities would be removed by the short time of travelling to this new end of the 
wire. Assuming below that one impurity creates a weak link between two parts of the wire 
and neglecting effects of other impurities one can apply the tunneling hamiltonian approach 
to describe the transport. In this approach current flowing through the weak link located, 
say, at x = inside the wire is given by the operator J(t) = — i[Tip^(0, t)ipp(0, t) — h.c), 
(e,h = 1), where T is the tunneling amplitude and ip^^x^) are the electron annihilation 
operators in the right < x < Lr and in the left —Ll < x < part of the wire, respectively. 



Then the average current under voltage V applied to the left lead is equal in the lowest 
order in T to: < J >= 2ii\T\ 2 J de[f(e — V) — /(e)]/?ij(e)pt(e — V), where / is Fermi 
distribution. The tunneling density of the right (left) end of the junction Pril){£) being 
the sum of the particle and hole densities can be extracted from the particle correlator as 
(1 - /(e))p(e) = 1/(2tt) Jdte iet < ^(0, t)^ + (0, 0) > or from the hole one as /(e)p(e) = 
1/(2tt) / dte iet < i}) + {0, 0)V>(0, t) >. 

To calculate the tunneling density to the right from the weak link pn let me first consider 
spinless fermions and apply bosonization to the ip field under condition of an elastic reflection 
from the boundary located at x — 0. ( Carrying out this calculations I will omit index "R" 



below. ) Above boundary condition known as "fixed" |14] in conformal field theory was 
sometimes addressed as "open" |15j in other considerations. Bosonic repersentation of the 
ip field reads i/j(x,t) = Y^ a =r,i i>a(%, t) = (27ra) _1 J2± exp{i(9(x,t) ± <f>(x, t))/2}, where ip r (i) 
is the right (left) going chiral component of if) and the 9 and fields are bosonic and 
mutually conjugated [9(x,t),<fi(y,t)] = 2nisgn(x — y). The elastic reflection means that 
i/ji(0,t) = e lS ip r (0,t) with an appropiate phase shift 5. This results in both: 

0(0, t) = 5, ^-d x 6(x, t) \ x=0 = W + (0, t)W(0, t) - ^ + (0, *)Vi(0, t) = 0. (1) 
Then the density of particle states could be found from 

pJe) = 1/(2tt) / dte iet < ^(0, t)^ + (0, 0) >= [ + °° dte *t+ll<o(o,t)e m >-<e* m >] (2) 

where the value of the free electron tuneling density was introduced as: po = (1 + 

COS(2<5))/(7T«). 

The problem reduces to finding the 9 field correlator. It can be done for the finite length 
piece of the wire adiabatically connected to the lead making use of the ITTL model [p|-|TT|j. In 
this model the Tomonaga-Luttinger interaction (J2 r ,iPa) 2 is switched on in the wire x < Lr 
and switched off outside. Then the Hamiltonian takes a bosonized form 



where function u(x) ensuing from the interaction can be approximated in the low energy 
limit by a step-function: u(x) — 1 if x > Lr and u(x) — u — 1/g < 1, otherwise. The 
correlator of the 9 field ordered in imaginary time T(x,y,r) =< T T 6(x,r)6(y,0) > can be 
shown to satisfy the following equation 

1 An 

{-mn% + 9 *} T & y> r ) = — - ^( r ) ( 4 ) 

V A U Z {X) V 

under the boundary conditions d x T(x,y,r)\ x= o = following from (|l|). Fourier transform 
of this correlator T(x, y, uj) is symmetrical under uj — > — uj. It can be compiled from the 
solutions of the homogeneous equation corresponding to Eq.(f|) 

2 2 a 

{ 2 2f n - ^}/«(^) = { 2 27 x - ^l^O) = — ^ - y) (5) 
An 

T(x,y,u) = [9(x - y)f u \x)K\m + 6(y - x) f u {y)h u [x)\ 

vW (uj) 

if these solutions meet boundary conditions: h'^O) = 0, fw(x) = exp(—ujx/v) at x — ► oo 
and positive o>. The Wronskian W(u;) is equal to — /^(0)/i w (0) and, hence, T(0, 0, cu) = 
— — l/(Z7i/ w (0))'. The only solution I need can be written as right going plus reflected left 
going waves at x < L R . The reflection amplitude r = — e~ 2v , tanh(?7) — 1/u equals minus 
reflection amplitude of the <fi field ||IT| , |T6|| due to duality symmetry and is negative for the 
repulsive interaction. Substituting this solution one can find T(0, 0, uj) = tanh(o;t LR + 77) 
with equal to the time of travelling from the junction to the right lead. Analytical 
continuation of this function [— T(0, 0, — iuj + 0)] brings us the value of the retarded Green 
function for the 9 field. Imaginary part of the latter multiplied by the Bose distribution 
function for holes 1 + /b(o;) and by a factor (—2) coincides with the Fourier transform of 
the correlator at uj what allows me to rewrite expression for the particle density of states 
(§) in dimensionless units as 

f\ P° f + °° a j • , r a , ; / Hmtan{uj + i V ) 
Pp {e) = ^J_ oo dpexp^ep + ul due \(1 + } (6) 

where the temperature and energy were scaled as e = t^e. The dimensionless cut-off 
parameter 7 equals [Ef^lr)^ 1 ■ The hole density of states Ph(e) can be found as Ph(e) = 



p p (—e). Below I will examine expression (|6]) at zero temperature when the expected effect 
of the interference on the density of states should be the most profound. The bosonic 
distribution factor 1 + fg restricts the integral to positive uj and the density of states at 
positive e coincides with the density of particle states. 
The form of the correlator used in (^) is equivalent to 

lt 2l l (7) 
7 + 2m + p j 

where p is dimensionless time p = t/tLR- This expression can be easily comprehended as 
if it is a product of the contributions of the 2n length paths connecting (0,p) and (0,0) 
points and undergoing n reflections from a x = Lr non-elastic boundary with the negative 
reflection amplitude r = — e~ 2r> and n reflections from the x = elastic boundary with 
unit reflection amplitude. Another form of the exponent in the left hand side of Eq. (||) 
exp{u(S(ip + 7) + / °° dz / (2z)[e~ zp — l][tanh(z +77) — tanh(z — 77)])} can be obtained after 
rotation of the integration contour uj = —iz. Here S(ip + 7) = —in f£ dse~^ v+l7T ^ 2 ^ s ~' l ' y,> (1 — 
e -i7r(s-i7)^-i ensues from the pole contributions. For a small 7] it gives expression for the 
correlator 

(8) 

where the divergent at p = An part and the smooth longly decaying tail finally approaching 
const/p at pn ^> 1 are separated. Such a long time decay brings about the finite p(0) value 
predicted before 1303 . 

Besides making use of the above form of the correlator in (0) it is convenient to represent 
the density as a sum p{e)/po = - L ^— (sin(7ra)r(l — m)^" -1 + 2r(e)) of a scale invariant 
density of the infinite long wire and the finite length correction. Then calculation of the 
latter reduces to Fourier transform of the integrable function at least if u < 3. It has been 
done numerically and results are depicted in Fig. la. They show that as the interaction 
is strengthened the dimensionless finite length correction to the density of states changes 



its period of oscillations from %/tLR for u = 1.4 to ir/(2tLR) for u = 1.86. At larger u the 
oscillations become more profound and practically cease to decay with increase of the energy. 
As u decreases to 1 the finite length correction r(e) goes to zero everywhere except for e = 
where the value r(0) approaches ir/2. Albeit the amplitude of the r function oscillations 
becomes much smaller than 1 in this limit it decreases slowly with increase of the energy. 

To generalize expression (H) for the density of states to the case of the spin electrons 
one can notice that the 9 a arising in Eq. © in this case may be represented as a sum over 
the charge and spin fields 6 a = (9 C ± 9 S ) / V2. Dynamics of the charge field is described 
by the Hamiltonian (H) as before while dynamics of the spin field is not affected by the 
interaction. It means that the spin generalization of the density of states requires change 
of im[tan(c<j + it])]/ tanh(?7) in (|^) into (1 + Jm[tan(u; + if])}/ tanh(r]))/2. The results of 
calculations for the spin electrons are shown in Fig. lb. They behave similar to the results 
of the spinless case, however, the interference structure is weaker than in the spin case at 
the same value of the interaction constant. 

Calculating the differential conductivity one can use the same form where the finite 
length correction is separated from the scale invariant infinite length contribution: d J / dV = 
^o 1 7 2(n_1) [(2w - l) v / vr2 1 - 2 "(r(M)r(l/2 + m))" 1 ^"- 1 ) + dj(v)}. Here Rq 1 is a free electron 
conductance of the junction and both the cut-off 7 and the applied voltage V are measured in 
the unit of the full traversal time = t lr^^ll so that v = Vt L . The finite length correction 
dj to the conductivity depends on relation t LR /t LL between the travelings times to the right 
and left leads. Its behavior is illustrated by Fig.2a,b. They show that as the voltage increases 
the correction grows up if u is larger than u ~ 2 and just oscillates otherwise. It is reasonable 
since the leading contribution to the correction comes from convolution of the oscillating 
r(e) function and the infinite length e u ~ l function. The amplitude of the oscillations behaves 
similar to the one of the density of states. It is on the order of the conductance for large u 
and smaller than conductance for u close to 1. However, in both cases the oscillations decay 
very slowly as the voltage increases what hopefully makes them observable. 

In summary, I have shown that the ITLL model accounting for the finite length of the 

6 



wire and predicting a finite zero temperature conductance also brings about the oscillating 
interference structure in the differential conductivity. The latter survives at energies much 
higher than the one corresponding to the wire length 1/if,. The period of the oscillations 
of the tunneling density of states corresponds to the right (left) traversal time for a weak 
interaction u ~ 1 and becomes two times shorter at larger u. The differential conductivity 
in turn oscillates with both periods of the right and left tunneling density of states. 

The author acknowledges N. Nagaosa, S. Tarucha for useful discussions. It is a special 
pleasure to thank T. Iitaka for his help in conducting calculations. This work was supported 
by the Center of Excellence and partially by the fund for the development of collaboration 
between the former Soviet Union and Japan at the JSPS. 
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FIGURES 



Dependences of the finite length correction (E F t L ) 1 ~ u 2r(E) /it to the density of states 
on the energy E measured in n over traversal time ti unit . In the spinless case 
(a) solid, dot dashed and dashed lines correspond to constant of interaction u equal to 
u = 1.396, 1.862, 2.825, appropriately. In the spin case (b) solid and dashed lines correspond 
to u = 1.862, 2.825. Long decay of the oscilations at u=1.396 is shown in the inset. 

Dependences of the finite length correction (EptL) 2 ^ 1 '^ 'dj(v) / Ro to the differential con- 
ductivity on volage v measured in n over traversal time ti unit. For spin electrons (a) 
u = 1.396 the solid line relates to the symmetrical position of the weak link, the dotted line 
to the case when relation between the lengths of the right and left shoulders is 1/4. In the 
spinless case (b) the lines 1,2 correspond to u = 2.164, 1.862, respectively. 
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